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I . I NTRODUCTION 



A jet, regardless of its latitudinal position in the at- 
mosphere, may contain regions of large vorticity gradients 
where the necessary condition for barotropic instability is 
locally satisfied. Synoptic-scale moving disturbances occur- 
ring at the level of a moderately strong easterly jet south 
of the Tibetan high near 200 mb level appear to arise from 
barotropic instability of the mean flow ( Kr i s hnamu rti , 1971a, 
1971b). If the observed disturbances are the result of baro- 
tropic instability, they will extract energy from the mean flow. 

Many investigators have studied the linear stability of 
barotropic zonal flows over the years. Kri shnamurti 1 s obser- 
vational and Colton's (1973) numerical studies both showed 
that the zonal variation of the jet apparently had significant 
effects on the dynamic behavior of the transient disturbances. 

Tupaz et al . (1978) further extended the study of the 
zonal variation of the basic flow. A numerical model was 
developed based on the linearized non-divergent barotropic 
vorticity equation on a beta-plane. This study examined the 
dynamics of transient barotropic waves in a region of variable 
mean wind. The mean zonal wind was an easterly hyperbolic 
secant-squared, or Bickley, jet and the mean meridional wind 
was derived so that the mean flow would be non-divergent. 

Using the numerical model, a periodic forcing was applied on 
the inflow (eastern) boundary to generate waves which were 
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allowed to propagate through the channel and out through the 
western boundary. As the waves moved, they reacted to the 
local stability properties of the mean flow, whereas at each 
point the fields varied pe ri o d i cal 1 y . The results, which 
were obtained from long term integration, were compared with 
the parallel flow theory by constructing a simple mechanistic 
analytical model which incorporated the local stability con- 
cept of the parallel flow theory. 

In this thesis, the model developed by Tupaz et al . (1978) 
is modified by changing the boundary conditions. The quanti- 
ties on the eastern and western ends of the channel are set 
equal to enforce cyclic continuity. This allows a wider 
variety of wind profiles to be considered as long as they 
are periodic in x. The use of these boundary conditions will 
also improve our understanding of the mathematical behavior 
of waves in a mean flow with downstream variation. 
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II . BASIC EQUATIONS 



A. GENERAL FORMULATION 

This study employs the barotropic vorticity equation to 
describe motion in a periodic channel. The development which 
follows Tupaz et al. (1978) begins with the barotropic vor- 
ticity equation on a beta (8) plane: 

H +u K +v K +6v , Q . DfC , (2.1, 



where 



_ 9 v 9 u 
5 ■ lx “ dy ' 



( 2 . 2 ) 



Here, 8 is the constant value of the north-south gradient of 
Coriolis parameter computed at 10° latitude. Q is a forcing 
function which represents the non-barotropi c and diabatic 
effects which are required to maintain the basic state vor- 
ticity field, and D* is a frictional coefficient. The mo- 
tion is assumed to be non-divergent so that a stream function 
(ijj) can be introduced which is defined by 



u = 



9 ^ 
9y 5 




The vorticity becomes 

? = V 2 ip . 



(2.3) 



(2.4) 
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The non-divergent barotropic vor ti ci ty equation may now be 
wri tten , 



9_ 

at 



V Z \p 



n nh + ii nh +g 11 = 

9y 9x 9x 9y 9x 



Q - D f V 2 ip . 



(2.5) 



The linearization of eq. (2.5) is accomplished by separat- 
ing the stream function into its mean (bar) and perturbation 
(prime) quantities of the form: 



<p( x ,y ,t) = (p(x,y) + f (x,y,t) , (2.6) 

where |^'| < < |p| . Substituting (2.6) into (2.5) and 
dropping the quadratic terms in ip' leads to 



a_ 

at 



|i vV 

9y 9x v 




dip 


" 17 


1— V 2 ip 
9x 9y v 


+ s &' - - 


D f v 2 ^ 



9 x 



V 2 p 



(2.7) 



It is assumed that Q' = 0 and Q is defined in such a way 
that jp is a steady-state solution of (2.5). Equation (2.7) 
can also be rewritten in Jacobian form such that 



v 2 II' = - jr (<mV ) - j (r ,v 2 *) -6 n' 

- D f vV , (2.8) 



where ( a , b ) 



9a 9b 9a 9b 
¥7 Ty ~ "9 y ~dx 



Equation (2.8) is a Poisson equation for the tendency of the 
perturbation stream function (-|r )• 



Equation (2.7) is solved in a channel defined by 



- D <_ y £ D 
X W 1 x 1 X E • 

The boundary conditions of the north-south walls y = +_ D 
are 

ip' = 0 . (2.9) 

Cyclic boundary conditions are used in the east-west direc- 
tion. Tupaz et al. (1978) applied a time-periodic forcing 
on the eastern boundary to introduce waves into the region 
from the east, and used a radiation condition for both ten- 
dency and vorticity at the western boundary to allow the 
waves to move smoothly through the boundary. 

B. ZONALLY VARYING BASIC FLOW 

The basic velocity field is defined by an easterly Bick- 
ley jet 

u(x,y) = U ( x ) sech2 (^f7y) " u 0 = ' |y * (2.10) 

Here, d(x) is a characteristic length scale of the jet, U(x) 
is the maximum velocity at y = 0, and U 0 is a constant 
velocity. The basic flow stream function is specified to be 
constant at the northern and southern boundaries (y = + D). 
Equation (2.10) can be integrated to give 

^ ( x , y ) = U ( x ) d(x){tanh( ^r^y ) + tanh ( d ( x j ) > 

- U 0 (y + D) + ij>(-D) , (2.11 ) 
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where 



$(D) + 2U D - iH-D) n 

U <*) " 2dfxJ coth * (2 * 12) 



Therefore, if d(x) varies slowly in x , so does the basic 
flow. The x-variation for the characteristic length scale is 
given by 

( X-X ) 

850 km + 350 km{cos [ 2 tt .. 0 ]} , 
d = L 

1200 km 



- x < x < x , 

o - - o 



x > x Q . (2.13) 



Here, is the longitude where the x-variation of the cosine 

function starts and L is the wavelength of this variation. 



III. FINITE DIFFERENCE EQUATIONS 



Equation (2.7) is solved by finite differences. Centered 
time differences are used except for the friction term which, 
is approximated with a forward step. The Jacobians are com- 
puted with the form which was developed by Arakawa (1966) for 
non-linear conservation. 

The tendency of the disturbance stream function is de- 
fined by 



T* . 
i > J 



= I^V 

' l 3t J i ,j 



(3.1) 



where superscript t is the current time and subscripts i 
and j refer to the x and y grid points, respectively. 
In the following finite difference equations, At(l hour) is 
the time step. Ax(375 km) and Ay(125 km) are the x and 
y grid point intervals, respectively. Equation (2.7) is 
written using the Leapfrog finite difference scheme: 

v2t ?,j * - -ji.jir.v 2 *) 



- 6 






I t 

*i-l 



2Ax 



-J- - D ! t- ! At 
f ,] 



i » j 



. , t- At , T . . 

+T i,j At 



(3.2) 

(3.3) 
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0h © hD @ I s ) 



The finite difference form for the Laplacian of the tendency 
is given by 



-2T^+T i _ , . , 

v 2jt s i +1 , .1 i , ,i i-1 . j + i , j+1 i ,.i i ,.i-1 



i > J 



Ax 



T. j ,,-2T. . +T 

' J ' 



and for the Jacobian terms it is given by 



(3.4) 



t 

i » J 



1 



4AxAy 



t( ^i+l , j ~ ^ i - 1 ,j )(? i!j + l' C i ,j-l ) 
^i , j+1 "^i , j-1 ^ i +1 , j~ ? i -1 , j ^ 



(3.5) 






i,j 4AxAy l + i +1 , j s i +1 , j+1 i +1 , j-1 



_ r 1 ^ ) 

’1 J.1 -l _ 1 / 



(3.6) 



^i-1 » j^i-1 , j+1 ~ ^ i - 1 , j - 1 ^ 






. t 



i , j+1 s i +1 , j+1 h i -1 , j+1 



. r.t ) 



+ ^ t ^ Uj 



. t 



i » j-1 i +1 , j-1 i-1 , j-1 



- i! 1 , ,_,)J . 



t 

i » j 



+ C 



i » j+1 ( ^i+l , j+1 


^i-1 , j + 1 


i .j-1 ( ^i+l .j-1 


' *1-1 .j-1 


i+1 ,j ( ^i+l .j+1 


*i+l .j-1 


i-1 , j ^i-1 , j + 1 


" ^i-l , j-1 
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(3.7) 



where 



® © O 

J 4 W.vVl -J* + J t + . (3.8) 

i » J 1 » J LJ j > -i 

3 

t 2 - 

Ji. . [^' ,V is expanded in the same manner as for (3.8). 

^ • J 

This Jacobian conserves both kinetic energy and enstrophy in 
the non-linear formulations (Arakawa, 1966). 

A forward time step is used for the initial time step and 
every 30 time steps thereafter to damp the internal gravity 
waves produced by the leapfrog method. A Matsuno finite dif- 
ference scheme is tested in place of the forward time step 
and the results show that the Matsuno scheme provides further 
damping, but that the change is not significant. 

The boundary conditions employed by Tupaz et al. (1978) 
relieved the problem of reflection of incident waves from the 
boundary back into the interior region (Pearson, 1974). An 
Euler-backward finite difference scheme was employed every 
time step to provide additional damping of the short waves. 

The frictional terms are evaluated at the previous time 
step in order to insure linear computation stability (Haltiner, 
1 971 ) . 

The Poisson Eq. (3.1) is solved for the tendency using the 
P0I5S subroutine which is a cyclic reduction algorithm for 
solving block tridiagonal systems of arbitrary dimensions with 
Dirichlet boundary conditions (Swarztrauber and Sweet, 1973). 

The subroutine POISS solves the linear system of equations 



of the form: 



Y . . 

1 > J 



A 1 X i-l , j +B i X i ,j tC i X i + l ,j +X i ,0-r 2X 1 , j +X i ,j+l 



for I = 1,2, ... M and J = 1,2, ... N . 

(3.9) 



The east-west boundary conditions for all quantities are 
represented by 



X 

o 

II 

X 

—t. 

w 

*— H 


i =0, I 




i , J + l = X i ,1 


+ 

* 

O 

II 

•*-> 


(3.10) 



It should be noted that Eq. (3.9) is not in the correct form 
which satisfies the model. Therefore, the model must be 
flipped prior to entering the solver and then flipped back 
upon exiting v/ith the solved eigenvalue-solutions. 
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IV. LOCAL STABILITY OF BASIC FLOW 



A parallel flow model of Williams et al. (1971) is used 
to determine the local growth rate of the variable mean flow. 
This is taken from Tupaz et al. (1978) and then transposed 
to fit the basic flow model used for this study in order to 
gain some insight on the stability character!' sti cs of the 
mean flow. This model is hereafter referenced as the par- 
allel flow model. If dip/dx is set equal to zero in (2.5) 
the governing equation of this model is 

‘If * 5 lx + “f 1 + 18 - H' " 0 ■ ' f 4 - 1 ) 

ay 

u(y) = -U sech 2 (^) + U Q 

where U is a specified constant that scales the magnitude 
of the central velocity of the Bickley jet (y=0). The 
characteristic length d is also a specified constant. 
Equation (4.1) is solved with the initial value technique 
used by Williams et al. (1971). Assuming that all perturba- 
tions are periodic in x, Eq. (4.1) is Fourier transformed in 
x and the resulting equations are integrated until the solu- 
tion grows. with a constant exponential rate. This approach 
gives the growth rate, phase speed and wave number structure 
of the most unstable meridional mode as a function of the 
zonal wave number k . 



20 



In general, these equations have a set of normal mode 
solutions and a continuous spectrum of solutions (Case, 1960; 
Pedlosky, 1964 and Yanai and Nitta, 1968). However, only 
the normal mode solutions can give significant growth and 
the most unstable mode will dominate after a sufficient 
period of time. 

Figure 1 shows the growth rate (n) corresponding to the 
most unstable wavelengths (L) as a function of x based on 
the parallel flow model. The parameters for the mean flow 
are given in Section V. The largest growth rate occurs at 
x=0 where the jet achieves its maximum central velocity 
[u = - 30 m sec - '*]. In the inflow and outflow region for 
x < - 11,100 km and x > 11,100 km, respectively, the growth 
rate becomes negative and approaches the linear frictional 
damping rate. The most unstable wavelengths range from 
3650 km at x=0 to 4600 km in the inflow and outflow regions. 

This variation follows the variation in the y-scale of jet 

d(x) which is given by (2.16) (Tupaz et al., 1978). 

Ei gensol uti ons of the most unstable discrete mode for 
two jet profiles are shown in Fig. 2. Only the lower half 
of the y domain is shown because the solutions are symmetric 
about y=0. Both solutions show barotropic growth as the 
phase tilt is opposite to the shear in u. The maximum tilt 
occurs at x=0 where the growth rate is the largest. In 
addition, the ei gensol uti on amplitude is a maximum aty=0 for 
x = 0, but at x = 9314 km the maximum has shifted to y = _+ 800 
km . 
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L (1000 km) 
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X (1000 km) 

Fig. 1. (a) The most unstable wavelength L(x), and (b) corresponding local growth 

rate n (x), based on the parallel flow model. 
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Parallel flow model wave structure of the most unstable mode for longitudes: 
(a) x=0 and (b) x = + 9619. 



V. NUMERICAL RESULTS FOR EXPERIMENT I 



The numerical investigation of the variable mean wind 
model comprises two experiments. Experiment I is the princi- 
pal one and it employs the same friction that was used by 

Tupaz et al. (1978). A friction coefficient of = 0.15 

- 5 - 1 

x 10 sec is used which is equivalent to an e- folding 
decay time of a pp roxi ma tel y 7.7 days. The second experiment 
is the frictionless case and will be discussed in Section VI. 
The following parameters are used in Experiment 1. 

2D = 4000 km, X E - x w = 40,1 25 km, u ( 0 , 0 ) = -30 m s -1 , 

U q = 0, L = 36,000 km. Ax = 375 km. Ay = 125 km , 

At = 1 hr, = + 18,000 km 

Figures 3 and 4 show the stream function ^(x,y) and the 
zonal velocity u(x,y) which are derived from the above 
parameters. The model has 108 grid intervals in the x-direc- 
tion and 32 intervals in the y-direction. Tupaz et al . (1978) 
used the same domain, but his mean flow was not the same at 
x =X E and x =x y However in this study, the scale of 
variation of the jet is reduced to insure that the basic 
zonal velocity would be periodic in the inflow and outflow 
regions. Equation (2.13) specifies that the jet scale, d(x), 
varies between 500 km at u ( x , 0 ) max and krn at ^ x ’ 0 ^ mi n ’ 

The initial disturbance field is of the form: 
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-21 -18 -15 -12 -9 -6 -3 0 

X (1000 km) -* 

Fig. 3. Experiment I. The ij»(x,y) field (x 10 7 m 2 sec 
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X (1000 kin) -* 

Fig. 4. Experiment I. The u(x,y) field (x 10 m sec 



(5.1 ) 



<!>' (x,y,o) 



sinTT^ cos 2tt IOx 
2D x £ -x w 



A wave number of 10 is chosen as a result of the average wave- 
length which was obtained by Tupaz et al . (1978). The initial 
conditions and the jet flow are both symmetric about y=0 and 
this symmetry is maintained throughout the integration. 
Therefore, only the lower half of the y plane will be depicted 
in the results of both experiments. 

The integration of the numerical model is extended for a 
period of 170 days. A close examination of wave structure, 
phase tilt, period, wavelength, and spatial growth is made 
for the segments 80-87 days, 125-133 days and 150-157 days. 
Overall growth occurs during the entire period and is expected 
since the waves travel through the region and then reappear 
at the input boundary. It was expected that the overall 
growth would be exponential in behavior. Figure 5 shows the 
root mean square for the entire domain of ip' plotted on a 
log scale vs time. If the growth were truly exponential, the 
curve in Fig. 5 would be linear. But in fact the curve is 
quasi-linear with departures from the linear behavior occur- 
ring periodically about every 56 days. Figure 5 shows the 
first departure occurs near t = 2 1 days and continues for 28 
days. It should also be noted that the growth rate is nega- 
tive from t=38 days to t=47 days. This departure from the 
linear appears again at t=75 days for 25 days and again at 
t = 1 3 5 days for 26 days. 
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ave 




Time (days) -*• 

Fig. 5. Experiment I. I^l ave as a function of time, where 
|ij>| a is the root mean square of ip'. 
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The i|>' field at t = 30 days is shown in Fig. 6. An entire 

train of baro tropic waves actually exists throughout the length 

of the channel, but the waves upstream of x = -3500 km and in 

the outflow region are not s hown . becaus e of their relatively 

small amplitude. This is the same type of wave behavior ob- 

\ 

served by Tupaz et al. (1978). 

The phase is calculated for the segments mentioned earlier. 
This is accomplished by observing each grid point in the ip ' 
field and recording the first three times it passes through 
zero in either direction. A seven day period (P) is used in 
order to insure that a full period is recorded. Tupaz et al. 
(1978) observed a periodicity of 3.25 days. The actual phase 
for each grid point is computed using the following formula: 

360° ( N- N ) 

0 - p — (5.2) 

where N represents the- first time it passed through zero and 
N q is the starting point in time of the sign calculations 
for each time segment. Although the wave continues to grow 
during the time segment, the period should not change appre- 
ciably. Figure 7 represents an analyzed section of the phase 
tilt for the segment 80-87 days. We note that all the waves 
shown tilt oppositely to the mean wind shear, which is neces- 
sary for barotropic instability. The maximum tilt is observed 

downstream of x=0 where the jet velocity is maximum. 

The periods at y = 0 and y =j^ 750 km for all three time seg- 
ments are shown in Figs. 8-10. Although the periodicity is 
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Fig. 6. Experiment I. The t|>' field at t = 30 days. 



270 180 90 360 270 u 180 u 90° 360 270 



O 




<- (un| OCCL) A 



31 



X (1000 km) -► 

Experiment I. Analyzed section of phase 6 during the segment 80-87 days. 
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Experiment I. Periods of the numerical model for the time segment 80-87 days 
for latitudes: (a) y = 0 and (b) y = + 750 km. 




Experiment I. Periods of the numerical model for the time segment 125-133 days 
for latitudes: (a) y = 0 and (b) y = + 750 km. 
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X (1000 km) 

Fig. 10. Experiment I. Periods of the numerical model for the time segment 150-157 days 
for latitudes: (a) y = 0 and (b) y = + 750 km. 



not constant, it is about 77 hours (3.2 days). This is also 
representative of the entire channel. Deviations near the 
inflow/outflow boundaries could be attributed to the small 
amplitudes in this region. Figures 10a and b show an increase 
in period at x = -6500 km to 93 hours (3.07 days) and 90 hours 
(3.75 days), res p ecti ve 1 y , and a decrease at x = 6500 km to 
39 hours (1.65 days) for a and b, respectively. Figure 8 
shows only an increase at x = 6500 km to 100 hours (4.17 days) 
and 93 hours (3.87 days) for a and b, respecti vely. This 
abnormal behavior could be connected to the departure of the 
growth rate from the linear behavior illustrated in Fig. 5. 

Figure 11 is an overview of the phase tilt for all three 
time segments. Although they are not the same, they are 
similar in the region from x = 3875 to x = -12000 km. They 
tend to be different in the region where the amplitude is 
small. 

The envelope of the wave packet is evaluated at y = -750 

km where the wave disturbance amplitude is large. The nor- 

malized \p' field is used because the waves are still grow- 
ing and results are something close to a sinusoidal variation. 
This envelope, <ip* (x)>, is obtained by recording the average 
of the magnitudes of the maximum and minimum V/|'H ave 
values that occur at each longitude over a period of 7 days. 

Here, l^l ave 1s the ro °t mean square ip' . Figure 12 shows 

the wave packet of all three time segments for y = _+ 750 km. 
Note that the maximum amplitude in Fig. 12 is near the point 
where the parallel flow growth rate drops to zero in Fig. 1. 
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X (1000 km) 

Fig. 11. Experiment I. Overview of phase tilt between y = 0 and y = + 750 km for the 

time segments: (a) 80-87 days (...), (b) 125-133 days ( ) and (c) 150-157 

days. 
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-21 -18 -15 -12 -9 -6 -3 0 3 5 9 12 15 

X (1000 kn) -* 

Fig. 12. Experiment I. Uave packet envelope, <ip' (x)>, of the numerical model for the 

time segments: (a) 80-87 days (...), (b) 125-133 days (.-•-) and (c) 150-157 
days ( — ). 



Although they are not the same, it should be noted that they 
are downstream of the max jet core. 

Tupaz et a 1 . (1978) derived the disturbance kinetic 
energy equation to better understand the dynamics of baro- 
tropic instability. They found a number of energy terms but 
observed that the Reynolds stress term, <(u'v' 3u/3y)>, was 
the only dominate source for disturbance kinetic energy. This 
supported the large amplitudes and spatial gradients of if/' 
encountered a considerable distance downstream of the jet. 

It is expected that this energy balance will also hold in this 
study. Since this is a region of tilt opposite to the shear, 
there is a conversion of mean flow kinetic energy to distur- 
bance flow kinetic energy. Figures 1 and 2 illustrate the 
large energy conversion of this type for the parallel flow 
model . 

Figures 13-16 show the phase angle 0 of the wave distur- 
bance as a function of y for three different time segments. 

The phase, 0 , for the most unstable wavelength of the par- 
allel flow model is also included (see Section 4). We note 
that the waves upstream of -10875 km tilt oppositely to the 
mean wind shear, which implies barotropic instability. In 
one case, the maximum tilt is downstream of x=0, which is 
where Tupaz et al . (1978) observed the maximum tilt. Between 
x = -10875 and x = 3750, the phase tilt compares closely to 
the parallel flow model 9 . At the eastern boundary (Fig. 15b), 
0£* compares rather well to 9 but 9-|* and 0^* show exces- 
sive tilt for the comparatively weak instability in this re- 
gion. This excessive fluctuation may partially be explained 



38 



(Phase) 



+ 



O 

CD 



O 



to 



I 



o 



LO 



o 

Osl 

I 




«- 000 L) A 



39 



Figure 13. Experiment I. Comparison of phase 9 of the numerical model for time 
segments (9^) 80-87 days (...), (e 2 *) 125-133 days (0 3 *) ISO- 

157 days ( — ) and the most unstable wavelength of the parallel flow 
model, 9 ( — ), for longitude x = 0. 
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Fig. 15. Experiment I. Comparison of phase 0 of the numerical model for time segments 
(0!*) 80-87 days (...), (0 *) 125-133 days (0 3 *) 150-157 days (— ) 

and the most unstable wavelength of the parallel flow model, 0 ( — ), for 

longitudes: (a) x = -10875 km and (b) x = 10875 km. 0 
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Fig. 16. Experiment I. Comparison of phase G of the numerical model for time 
segments (6-,*) 80-87 days (...), (e 2 *) 125-133 days (G 3 *) 15 

157 days ( — ) and the most unstable wavelength of the parallel flow 
model, 0 ( — ), for longitude x= -21000 km. 



by the small amplitudes which are present in this region. 

Near the outflow boundary (Fig. 16) the tilt is reversed in 
★ * 

and 0 ^ . This indicates dynamic stability or a flow of 
energy from the disturbance back to the mean flow. This can 
also be detected in Figs, lib and c. In almost all cases, we 
observe that the tilt of the wave disturbance qualitatively 
adjusts to the local stability of the mean flow. Since the 
parallel flow model can only solve the most unstable discrete 
mode, dynamic damping is not indicated in Fig. 15 for the 
outflow region. Note that the phase tilt is unstable on the 
average. This is because the waves show a net growth on 
moving through the region. 

the wavelengths for the numerical model are computed 
using the fo rmu 1 a : 



2tt 2Ax 
A0 



(5.3) 



where A0 represents the change in phase over 2Ax . 
Equation (5.3) can be simplified and expressed in degrees 
resulting in 



360° 2Ax 
A0° 



(5.4) 



In addition, a smoothing process is applied in those areas 
where the wavelength varies noticeably to give a more accurate 
picture of the wavelength throughout the channel. This cal- 
culation is made by choosing a questionable point ( O ) in 
the y domain of the phase field and adding 180° to it. This 
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new phase is then located by moving in both x directions. The 

distance between the new phases is calculated in terms of the 

number of x grid points and multiplied by Ax to give the 

corrected, or smoothed, wavelengths. Figures 17-19 show the 

★ 

variation of wavelength, L(x) , for all three time segments 
at y = _+ 750 km. L(x) q represents the smoothed wavelength. 
Comparing the wavelengths of y = 0 (not shown) and y = _+ 750 km, 
it is observed that the wavelengths are quite similar for 
each time period. The maximum wavelengths are within +_ 4500 
km of x =0 and tend to decrease toward the boundaries when 
the wavelength pattern becomes noisy. This is the same region 
that the phase tilt (Figs. 15b, 16) experiences large devia- 
tions from the parallel flow model due to the small ampli- 
tudes. Tupaz et al. (1978) observed that the wavelength at 
y = _+ 750 km was larger upstream and smaller downstream of 
x =0 than at y= 0 . 

★ 

Figures 20-22 contain the spatial growth, m.j _ ^ , computed 

from the numerical model, res pecti vel y . Figure 23 is the 

★ . 
smoothed version of m ^ _ 3 and the spatial growth rate (m) of 

the parallel flow model obtained by using 

m = -n/C r (5.4) 

where C r is the phase of the waves and n is the local 

★ 

growth rate (Tupaz et al., 1978). The value of m^ ^ 1 s com- 
puted directly from the wave packet envelope. In Fig. 22, 

★ 

we see that m^ ^ has a larger maximum than m and the maximum 
is shifted slightly upstream from the jet maximum. All three 
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Experiment I. Latitude y = + 750 kin (a) Wavelength L(x)£* for time 
segment 125-133 days, (b) Wavelength L(x) 9 computed by smoothing L(x) 
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X (1000 km) -* 

Fig. 19. Experiment I. Latitude y = + 750 km. (a) Wavelength L(x)^,* for time 
segment 150-157 days, (b) Wavelength computed by smoothing L(x) 
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Fig. 20. Experiment I. Spatial growth rate, for time segment 80-87 days. 
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Experiment I. Spatial growth rate, m(x) *, for time segment 125-133 days. 
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Fig. 22. Experiment I. Spatial growth rate, m(x)^*, for time segment 150-157 days. 
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curves pass through zero at the same point in which the wave 

packet envelopes are a maximum (see Fig. 12). Note again the 

noisy pattern near the eastern boundary. Tupaz et al . (1978) 

★ 

observed these same features except that the m curve was 
skewed to the left with respect to m and the jet maximum. 

This was the result of the tilt structure in the wave lagging 
spatially behind the expected value from the local stability 
conditions. The (*) represents the maximum and minimum ob- 
served by Tupaz et al. (1978). Note the similarity of the 

★ ★ 

maximum and minimum to and m-j respectively. The reason 

★ 

that m.j ^ shows more damping in the outflow region than the 
parallel flow model is because the latter can only give the 
eigensolutions of the most unstable mode. 
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VI . NUMERICAL RESULTS FOR EXPERIMENT II 



Experiment II is the frictionless case and employs 
the same parameters and calculations procedures that are 
used in Experiment I. However, although the integration for 
this experiment is extended for 170 days, only the time seg- 
ments of 80-87 days and 125-133 days are included in these 
resul ts . 

Figure 24 shows the root mean square ip' for the entire 
domain. Note that the total growth is nine orders of magni- 
tude greater than that observed in Experiment I (see Fig. 5). 
In addition, the curve in Fig. 24 is more linear in appear- 
ance than Fig. 5 but still experiences the departures from 
the linear behavior during the same time interval that were 
observed when friction was included. 

The periods at y=0 and y = + 750 km show very little 
variation and closely resembled those pictured in Figs. 7 
and 8. The predominating period which is representative for 
the entire channel is the same as Experiment I (77 hours or 
3.2 days ) . 

Figure 25 is an overview of the phase tilt for the two 
time segments. The region of similarity extends from x = 

5000 km to x = -15,000 km. Uhen the amplitudes are small, 
the behavior becomes erratic as observed in Experiment I 
(see Fig. 11). 
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Fig. 24. Experiment II. |i|>| g as a function of time, where M ave 
is the root mean square of i|>' . 
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X (.1000 km) + 

Fig. 25. Experiment II. Overview of the phase tilt between y = 0 and y = + 750 km 
for the time segments: (a) 80-87 days (...) and (b) 125-133 days T ). 



The wave packet envelope, <ty'(x)>, is shown in Fig. 26. 
When comparing Fig. 26 to Fig. 12, it is observed that the 
curves are almost identical. Since the tilt is opposite to 
the shear, there is a conversion of mean flow kinetic energy 
to disturbance flow kinetic energy. 

Figures 27-30 show the phase angle tilt of the wave dis- 
turbance y-structure compared to the most unstable wavelength 
of the parallel flow model 0 . In both cases, the maximum 
tilt was at x=0. Comparing them to Figs. 13-16 in Experiment 
I, we see a similar structure except for Figs. 27 and 29a 
which show greater tilt between 500-750 km. Note again the 
excessive tilt of 0^* for the comparatively weak instability 
shown in Fig. 29b, and the elimination of the growth of 6^* 
of y greater than -1000 km in Fig. 30 for the Figs. 15b and 
16, respectively, in Experiment I. The curve for 0^* stresses 
the dynamic stability which can also be observed in Fig. 25b. 

The wavelength patterns for y=0 and y = + 750 are very 
similar to those observed in Experiment I. The maximum wave- 
length occurred within + 4000 km of x=0. 

The spatial growth rates (Figs. 31-33) are near mirror 
images of Figs. 20, 21 and 23. The growth rates are a little 
larger than in Experiment I as would be expected because there 
is no friction. 
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Fig. 26. Experiment II. './ave packet envelope, <i/i'(x)>,of the numerical model for 
the time segments: (a) 80-87 days and (b) 125-133 days. 
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Fig. 27. Experiment II. Comparison of phase 0 of the numerical model for time 

segments (e^) 80-87 days ( . . . ) , (e-*) 125-133 days ( ) and the most 

unstable wavelength of the parallel flow model, 6 ( — ), for lonqitude 

v _ r\ f> ^ 



(Phase) + 0° (Phase) 



o 



l n 



o 





0001) A 



59 



Fig, 28. Experiment II. Comparison of phase 0 of the numerical model for time 

segments (0-, ) 80-87 days (...), (0 2 ) 125-133 days ( ) and the most 

unstable wavelength of the parallel flow model, 0 ( — ), for longitudes: 
(a) x = -3750 km and (b) x = 3750 km. 0 
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Experiment II. Comparison of phase 0 of the numerical model for time 
segments (0-j*) 80-87 days (...), (0p*) 125-133 days (•-•-) and the most 
unstable wavelength of the parallerflow model, 0 ( — ), for longitudes: 

(a) x = -10,875 km and (b) x = 10,875 km. u 
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Fig. 30. Experiment II. Comparison of phase 0 of the numerical model for time 
segments ( 0 -|*) 80-87 days (...), ( 0 2 *) 125-133 days (<-•) and the most 
unstable wavelength of the parallel flow model, 0 ( — ), for longitude 
x = -21 ,000 km. 0 
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Fig. 31. Experiment II. Spatial growth rate, m(x) 1 *, for time segment 80-87 days. 
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Fig. 32. Experiment II. Spatial growth rate, m(x)p*, for time segment 125-133 
days. 
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VII. SUMMARY AND CONCLUSIONS 



This study solved the model developed by Tupaz et al. 
(1978) with different boundary conditions. They examined 
the behavior of linear waves in a zonally va ry ing easterly 
jet which was barotropically unstable in certain longitudes. 
The waves were forced with a period w on the eastern 
boundary and a radiation condition was applied on the west- 
ern boundary to allow the waves to move smoothly out of the 
region. After a period of adjustment the solution at every 
point varied with period u) , but as each wave moved through 
the domain it experienced growth or decay in response to 
the local mean fTow. 

In the present study, the Tupaz model was modified by 
changing the boundary conditions to cyclic continuity on the 
eastern and western ends of the channel. In addition, the 
mean was modified slightly so that it satisfied cyclic con- 
tinuity. Experiment I corresponded to Experiment I in 
Tupaz et al. (1978) and Experiment II was frictionless, but 
otherwise the same as Experiment I. 

In Experiment I, the root mean square stream function 
grew with time because as each wave moved through the channel, 
its net growth was added when it reentered the region. This 
growth was appro xi mate! y exponential with the departures from 
exponential growth having a period of about 56 days. In order 
to compare the wave behavior in this study with that of 
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Tupaz et al . ( 1 978 ), the fields were normalized to remove 
the overgrowth in time. A procedure was developed to deter- 
mine the phase and amplitude of the waves during a particular 
time segment. The following time segments were chosen for 
analysis: 80-87 days, 125-133 days, 150-157 days. The wave 

structure features were not the same for the three periods, 
but the features were very similar in the large amplitude 
regions. The predominant period was about 3.2 days which was 
very close to the forced period of 3.25 days used by Tupaz 
et al . (1978). The largest phase tilt was found near the 
jet maximum and in two of the segments the phase tilt reversed 
in the stable regions as was found by Tupaz et al . (1978). 

The amplitude of the wave envelope was a maximum near the 
point where the parallel flow growth rate was zero. The 
average of the maximum spatial growth rates was very close 
to the value obtained by Tupaz et al. (1978). This verified 
their result that the spatial variation in the growth rate 
could augment the maximum growth rate. The downstream wave- 
length was somewhat erratic, but it was generally a maximum 
near the jet maximum and it was a minimum where the jet was 
weak. The results for the frictionless case were similar, 
but the total growth was much larger. 

In general, the wave behavior with cyclic boundary condi- 
tions was very similar to that obtained by Tupaz et al . 

(1978). However, the solutions in this thesis were more dif- 
ficult to analyze because the wave structure in the three 
time segments was not the same. This was undoubtedly related 
to the non-exponential overall growths. Apparently as the 
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waves came out of the baro tropic damping region and entered 
the jet region they did not have the structure to begin growth 
smoothly. In this region the wave phase structure was par- 
ticularly erratic. When the present model was used to deter- 
mine wave structure and behavior with other wind profiles, 
the structure should have been computed frequently and then 
the typical structure could have been obtained by averaging. 

This study should be extended to the nonlinear case, by 
including a forcing term to maintain the mean flow when there 
are no waves present. When waves are introduced the modifi- 
cation of the mean flow could be studied. 
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